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I.  INTRODUCTION 


The  superior  energy  dissipating  efficiency  of  carbon  materials  in 
ablative  thermal  protection  system  applications  is  well-known.  However,  in 
advanced  systems  applications  such  as  planetary  entry,  both  the  solid-vapor 
interphase  mass  transport  rate  and  the  carbon  (graphite)  surface  temperature 
may  achieve  levels  considerably  higher  than  in  any  previous  use.  Thus,  ques¬ 
tions  concerning  energy  dissipation  efficiency  in  these  new  regimes  must  be 
addressed  if  survivable  systems  are  to  be  designed. 

Historically,  carbon  ablation  models  have  been  based  on  the  assumption 

1  2 

of  equilibrium  thermodynamics,  *  i.e.,  the  carbon  vapor  state  at  the  abla¬ 

tion  surface  is  assumed  to  be  given  by  equilibrium  thermodynamic  properties  at 

the  surface  temperature.  At  temperatures  above  40Q0°K,  the  equilibrium  total 

2  3 

vapor  pressure  of  carbon  is  presently  uncertain,  with  some  literature  values  * 

4 

being  as  much  as  an  order  of  magnitude  greater  than  the  JANAF  predicted  value. 
Also,  reported  melt  temperatures  range  from  about  4000  to  4400°K.^  When  melt¬ 
ing  occurs,  the  ablation  efficiency  may  decrease  by  as  much  as  a  factor  of 
three  if  the  liquid  carbon  is  stripped  away  without  vaporizing.  Thus,  even 
from  the  perspective  of  an  equilibrium  ablation  model,  carbon  ablation  effi¬ 
ciencies  above  4000 °K  cannot  at  the  present  time  be  predicted  with  desired 
accuracy. 


^Dolton,  T.  A.,  H.  E.  Goldstein,  and  R.  E.  Mauer,  "Thermodynamic  Performance 
of  Carbon  in  Hyperthermal  Environments,"  Progress  in  Astronautics  and  Aero¬ 
nautics:  Thermal  Design  Principles  of  Spacecraft  and  Entry  Bodies,  Voi.  71, 

Ed.  J.  T.  Bevans,  Academic  Press,  NY,  l & ,  pp  169-201. 

%aker,  R.  L.,  "Graphite  Sublimation  Chemistry  Nonequilibrium  Effects,"  AIAA 
Journal ,  Vol.  15,  Oct.  1977,  pp.  1391-1397.  - 

3Kratsch,  K.  M.,  et  al.,  "Graphite  Ablation  in  High  Pressure  Environments," 

AIAA  Paper  68-1153,  1968. 

^ JANAF  Thenmochemical  Tables,  National  Bureau  of  Standards,  NBS-37,  June  1971. 

5Baker,  R.  L.  and  P.  G.  Crowell,  "Graphite  Material  Ablation  Performance  in 
High  Thermal  Radiation  Environments,"  Progress  in  Astronautics  and  Aeronautics: 
Entry  Heating  and  Thermal  Protection,  voiT  by,  Ed.  w.  oistad,  AiAa,  NV,  ia8U, 


This  picture  is  further  complicated  by  the  multi-species  character  of 

C 

the  carbon  vapor  state  and  the  nonequilibrium  nature  of  carbon  ablation  in 

2  5  7 

advanced  applications.  »  *  To  encompass  these  features,  an  ablation  model 
must  express  the  interphase  mass  transport  rate  in  terms  of  the  vaporization 
kinetics  of  each  carbon  species. 

In  such  a  nonequilibrium  model,  these  rates  are  given  in  terms  of  the 

species  vaporization  coefficients.  Carbon  species  vaporization  coefficients 

have  been  experimentally  determined  only  at  the  relatively  low  temperature  of 
6  8 

2700*K.  *  There  is  apparently  no  published  information  concerning  their 
temperature  dependence. 

In  summary,  the  primary  thermochemical  properties  required  in  order 
to  accurately  predict  expected  carbon  mass  loss  rates  per  unit  input  of  energy 
are  the  vapor  pressure,  the  melt  temperature,  and  the  kinetic  rates  associated 
with  a  nonequilibrium  phase  change.  The  current  uncertainty  or  nonexistence 
of  values  for  these  properties  at  operating  temperatures  above  4000oK  pro¬ 
vided  the  motivation  for  the  work  reported  herein.  The  principal  objectives 
were  to  determine  the  above  properties  through  a  combined  experimental  and 
analytical  effort. 

The  vapor  pressure  and  vaporization  coefficient  determination  tech¬ 
niques  used  in  this  work  represent  a  significant  departure  from  standard  high 
temperature  methods,®’^  i.e.,  Knudsen  cells  and  free  vaporization  experi¬ 
ments.  The  present  method  is  based  upon  correlating  experimental  carbon  laser 
ablation  data  with  an  analytical  model  of  the  nonequilibrium  phase  change 


6Palmer,  H.  B.  and  M.  Shelef,  "Vaporization  of  Carbon,"  Chemistry  and  Physics 
of  Carbon,  Vol.  4,  Ed.  P.  L.  Walker,  Jr.,  Marcel  Dekker,  Inc.,  NY,  1968. 
'Baker,  R.  L.,  "An  Irreversible  Thermodynamics  Model  for  Graphite  Sublimation 
In  Radiation  Environments,"  Progress  in  Astronautics  and  Aeronautics:  Outer 
Planet  Heating  and  Thermal  Protection  Systems,  Vol.  64,  Ed.  R.  Viskanta,  AIAA. 

8V,  pp.  2\U-M7.  - - 

Zavltsanos,  P.  0.,  "The  Vaporization  of  Pyrolytic  Graphite,"  General  Elec¬ 
tric  Co.,  GE  R  66S031,  May  1966. 
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process.  The  biggest  advantage  of  this  new  approach  Is  that  It  allows  the 
upper  temperature  limit  of  the  carbon  data  to  be  extended  from  the  normal 
2700*K  up  to  4600*K.  A  unique  approach  for  the  determination  of  carbon  melt¬ 
ing,  based  on  detecting  changes  in  surface  reflectivity,  has  also  been  used. 

Experimental  carbon  ablation  data  have  been  obtained  over  a  wide 
range  of  surface  temperature  (3300  to  4600#K)  by  employing  a  continuous  wave 
laser  for  data  at  the  lower  end  of  this  temperature  range  and  a  pulsed  laser 
for  achieving  higher  temperatures.  A  nonlinear  analytical  model  describing 
the  nonequilibrium  gas  dynamics  and  chemical  rate  processes  that  occur  when 
carbon  is  rapidly  vaporized  has  also  been  developed10  as  an  extension  of 

pc  7 

earlier  work  with  linearized  models.  ,3».  This  report  discusses  comparison 
of  model  predictions  with  the  experimental  data  and  conclusions  from  this 
regarding  the  desired  high  temperature  thermochemical  properties  of  carbon. 

The  supporting  analytical  work  is  described  in  Section  II,  followed 
by  a  discussion  of  data  interpretation  in  Section  III.  Conclusions  are  sum¬ 
marized  in  Section  IV. 


^Margrave,  J.  t.,  Ed.  The  Characterization  of  High- Temperature  Vapors, 

John  Wiley  &  Sons,  NY,  '1%7. 

'°Baker,  R.  L.,  "Carbon  Nonequilibrium  Phase  Change,”  Office  of  Naval  Research 
Interim  Report,  TR-0081 ( 6728-02 ) - 1 ,  The  Aerospace  Corporation,  El  Segundo,  CA, 
December  1981. 
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II.  SUPPORTING  ANALYSES 


In  order  to  obtain  the  desired  thermochemical  properties  from  the 
experimental  carbon  laser  ablation  data,  an  analytical  description  of  the  non¬ 
equilibrium  gas  dynamics  and  chemical  kinetics  of  the  phase  change  process  is 
needed.  The  nonlinear  Knudsen  layer  model  which  we  have  developed  from  this 
purpose  has  been  described  in  Ref.  10.  In  Subsection  II-A  below,  the  model  is 
briefly  reviewed  and  recent  improvements  to  it  are  described.  Following  this, 
baseline  model  predictions  and  the  sensitivity  of  calculated  results  to  model 
parameters  are  discussed.  In  Subsection  II-B,  complementary  analytical  work 
in  progress  through  a  joint  effort  with  the  University  of  Texas  at  Austin  is 
descri bed. 

A.  NONLINEAR  KNUDSEN  LAYER  MODEL 


1.  Model  Description 

As  shown  schematically  in  Fig.  la,  the  Knudsen  layer  is  a  thin  region 
(the  order  of  a  few  mean  free  paths  in  thickness)  immediately  above  the  rapidly 
subliming  or  vaporizing  carbon  surface.  In  this  region,  nonequilibrium  gas 
dynamic  and  chemical  rate  processes  occur  which  determine  the  carbon  vapor 
state  as  it  is  transported  away  from  the  ablating  surface.  Macroscopic  mea¬ 
surements  of  pressure  p,  density  P,  temperature  T,  and  convective  velocity  u 
can  only  be  made  in  the  continuum  gas  dynamics  region  above  the  Knudsen  layer 
edge.  A  Knudsen  layer  model  provides  the  relationships  between  these  measur¬ 
able  variables  and  the  equilibrium  (saturation,  subscript  s)  thermodynamic 
properties  of  the  carbon  vapor  at  the  solid  or  liquid  surface. 

The  modeling  approach  we  have  used  represents  an  extension  of  the 
method  first  described  by  Anisimov.^  It  is  based  upon  equations  for  the 


^Anisimov,  S.  I.,  "Vaporization  of  Metal  Absorbing  Laser  Radiation,"  Soviet 
Physics  (JETP) ,  Vol.  27,  No.  1,  1968,  pp.  182-183. 
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conservation  of  mass,  momentum,  and  energy  across  the  Knudsen  layer.  This  is 
shown  schematically  in  Fig,  lb,  where  J  may  be  thought  of  as  a  flux  of  mass, 
momentum,  or  energy.  At  the  edge  of  the  Knudsen  layer,  the  fluxes  are  repre¬ 
sented  by  their  well-known  expressions  from  continuum  gas  dynamics.  At  the 
surface,  ayJs  represents  directional  fluxes  away  from  the  surface  due  to 
a  change  of  phase  from  the  solid,  and  Jb  represents  incoming  directional 
fluxes  caused  by  back scattering  of  molecules  during  collision  processes  within 
the  Knudsen  layer.  A  fraction  of  Jb,  acJb,  is  condensed  while  the 
remaining  fraction  (1  -  <*c)Jb  is  reflected.  The  quantities  ay  and  ac 
are  the  vaporization  and  condensation  coefficients,  respectively. 

Following  Anisimov  and  Bird,  the  mass,  momentum  and  energy  conserva¬ 
tion  equations  across  the  Knudsen  layer  for  carbon  species  C^,  1*1 ,  2,  ...n, 
may  be  written 


VT 


n 

E 

i=l 


(2-  oc  )  41  f*« 

I  ci  Mi 


x<) 


(2) 


12B1rd,  G.  A.,  Molecular  Gas  Dynamics,  Oxford  University  Press,  London,  1976. 
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The  functions  f(x^)  are  given  by 

_x2 

f 1  { )  ■  rr^^x^erfc(x^)  -  e  *  (4a) 

X.  -X? 

Axj)  *  (X1  +  ?)erfc{xi)  ■  ~^T7Z  e  1  {4b) 

f3(Xi)  =  x^  (x2  +  1)7 r  1/2  erfcfx^ )  -  (x2  +  2)e  1  (4c) 


where 


(4d) 


quantities  not  yet  defined  in  the  above  equations  are  the  molecular  weight  M, 
the  convective  velocity  w  of  the  backscattered  molecules  striking  the  surface, 
and  G  which  is  related  to  the  ratio  of  specific  heats  8  by 


« 


The  quantities  In  Eqs.  (1)  through  (4)  are  dimensional  except  for  x^, 
Tb,  and  w.  The  latter  two  variables  are  nondlmenslonalized  by  T  and  u,  respec¬ 
tively,  and  Xj  is  directly  proportional  to  a  species  Mach  number  M^,  i.e. 


Tb  ■ 

VT 

(6a) 

w  = 

w/u 

(6b) 

xi  = 

Vv2  "iVi 

(6c) 

When  it  is  assumed  that  T{ 
equations  reduce  to  those 

&  =  T  and  w  =  u,  then  ?b  =  w  =  1  and  the  above 
described  in  Ref.  10.  The  three  conservation 

equa- 

tlons  do  not  allow  one  to  solve  for  Tb  and  w.  An  ad  hoc  assumption  must  be 
made  in  order  to  solve  for  the  other  quantities.  Anisimov  assumed  ?b  =  w  =  1. 
If  the  equations  are  written  as  above,  these  quantities  may  be  parametrically 
varied  over  their  entire  allowable  range,  i.e. 


1  <  TK  <  T  /T  ,  0<  w  <  1  (7) 

b  s 

In  this  way,  the  sensitivity  of  the  calculated  results  to  the  ad  hoc  closure 
assumption  can  be  assessed. 

2.  Baseline  Calculated  Results 


In  all  of  the  calculations  in  this  report,  it  has  been  assumed  that 

a  ■  a  *  a.  Thus,  from  this  point  on  the  subscripts  c  and  v  will  be 
C  V 

dropped.  Additional  assumptions  must  be  made  before  the  coupled  nonlinear 
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equations  (1)  through  (3)  can  be  solved.  These  assumptions  are  summarized 
below: 


Closure  assumption  for 
conservation  equations 


Baseline  Sensitivity  Study 

Tb  =  w  »  1  1  <  Tfa  <  0  <  w  <  1 


Knudsen  layer  chemical 
species  relaxation 


Frozen 


Equilbrium 


Knudsen  layer  internal 
energy  states  relaxation 

Species  vaporization 
coefficients 

(adsorbed  layer  kinetics) 


Equilibrium  Frozen 

Frozen  Equilibrium 


The  baseline  assumptions' are  supported  by  the  relatively  large  number 
of  molecular  collisions  (thousands)  required  for  species  relaxation,  the  small 
(5  to  10)  number  of  molecular  collisions  required  for  rotational  energy  states 
to  relax,  and  the  short  residence  time  in  the  surface  adsorbed  layer  for  mole¬ 
cules  passing  from  the  solid  to  the  continuum  vapor  when  the  surface  tempera- 

13 

ture  is  3000°K  or  greater.  When  the  kinetic  processes  in  the  adsorbed 
layer  are  frozen,  the  species  vaporization  coefficients  are  assumed  to  be 
given  by 


a,-  S  a 


(8a) 
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and  when  these  processes  are  in  equilibrium,  it  can  be  shown  ’  that 


«i  -  («.,) 


(8b) 


13Rosenblatt,  G.  M.,  Personal  Communication,  Pennsylvania  State  University, 
Department  of  Chemistry,  University  Park,  PA. 

'^Klelnschmldt,  P.  D.  and  G.  M.  Rosenblatt,  "Relative  Vaporization  Rates  of 
A$2(g)  and  AS4(g)  from  an  Arsenic  (111)  Surface,"  American  Society  of  Mass 
Spectrometry  Meeting,  June  1977. 
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Baseline  calculated  results  are  shown  in  Fig.  2.  The  ratio  of  quanti¬ 
ties  across  the  Knudsen  layer  (see  Fig.  la)  is  presented  as  a  function  of  the 
Mach  number  M  -  u/  x/» RT4K  at  the  Knudsen  layer  edge.  The  parameter  variation 
with  a  is  also  shown.  From  this  figure,  the  values  of  p/ps  for  a  Mach  number 
of  1,  i.e.,  p*/p$,  can  be  obtained.  The  ratio  of  stagnation  pressure  pQ  to 
ps  can  then  be  determined,  since  for  M  =  1 


Po 

P* 


<^f 


T 


Po 

Ps 


StL* 

Ps  P1 


(9) 


Thus,  the  desired  relationships  between  quantities  measured  in  the  laser  abla¬ 
tion  experiments,  i.e.,  the  mass  loss  rate  m  =  pu,  pQ,  and  T  ,  and  the 
vapor  pressure  p$  and  vaporization  coefficient  a  have  been  established. 

The  rather  weak  nonlinear  dependence  of  P0/Ps  and  m*/m  on  a  can  be  seen 
by  renormalizing  and  plotting  these  quantities  as  shown  in  Fig.  3.  Note  that 
m*  is  the  value  m  when  the  Knudsen  layer  edge  Mach  number  M  is  1  and  ms  = 
PSJKTS/2*M  .  Thus,  Fig.  3  represents  calculated  solutions  of  Eqs.  (1) 
through  (3)  as  a  function  of  a  for  the  case  M  =  1. 

I 

Predicted  results  such  as  those  shown  in  Fig.  3  can  be  used  directly 
to  obtain  vaporization  coefficients  from  the  experimental  data.  The  method  is 
described  further  in  Subsection  III-B-1.  Similarly,  predicted  results,  such 
as  those  presented  in  Fig.  2,  can  be  used  to  obtain  vapor  pressure  ps  from 
experimental  data  as  described  in  Subsection  III-B-2  and  Ref.  10. 

3.  Sensitivity  Studies 


a.  Closure  Assumption  for  Conservation  Equations 

The  sensitivity  of  calculated  results,  for  the  primary  quantities  of 
direct  interest  p/p$  and  m/m$,  to  the  closure  assumption  Tb  =  w  =  1  is  shown 
in  Fig.  4.  Both  Tb  and  w  were  systematically  changed  over  their  entire  range 
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Fig.  2.  Baseline  Calculated  Results 


Closure  Assumpti 


of  possible  values,  l.e.,  Eq.  (7).  As  noted,  p/p$  Is  extremely  Insensitive  to 
the  closure  assumption.  The  other  variable  m*/iiis  Increases  by  a  maximum  of  5 
percent  when  It  Is  assumed  that  w  *  0  and  ■  T$/T.  Physically,  assuming 
w  ■  0  Mould  not  appear  to  be  realistic  when  the  Hach  number  Is  large.  Thus, 
especially  for  M  a  1,  the  baseline  calculated  results  for  m*/ms  should  be  very 
close  to  the  correct  value. 

b.  Chemical  Species  Equilibrium  In  the  Knudsen  Layer 

If  the  carbon  species  concentrations  In  the  Knudsen  layer  are  assumed 
to  be  In  chemical  equilibrium  rather  than  frozen  at  their  saturation  condition 
values,  as  In  the  baseline  calculations,  the  predicted  results  change  as  shown 
In  Fig.  5.  The  primary  quantities  of  Interest,  p/p$  and  (m/ms)  as  well  as 
P/Ps,  are  relatively  Insensitive  to  the  rate  of  chemical  relaxation  processes 
In  the  Knudsen  layer.  The  temperature,  however.  Is  more  sensitive  to  this 
assumption,  since  equilibration  Involves  changes  In  Individual  species  concen¬ 
trations  with  relatively  large  net  energy  release. 

c.  Frozen  Internal  Energy  States  In  the  Knudsen  Layer 

The  baseline  calculations  assume  that  the  Internal  energy  states  of 
the  polyatomic  carbon  species  are  In  equilibrium  throughout  the  Knudsen  layer. 
If,  Instead  of  this,  the  energy  In  these  states  Is  frozen  at  the  saturation 
temperature  T$  value,  then  much  less  energy  is  available  to  be  transferred 
to  translational  and  directed  (convective)  velocity  states.  As  a  consequence, 
when  these  states  are  frozen  the  temperature  at  the  edge  of  the  Knudsen  layer 
Is  considerably  reduced. 

This  effect  and  the  effects  on  m,  p,  and  p  are  seen  in  Fig.  6.  For 
a  Knudsen  layer  edge  Mach  number  of  1,  both  the  temperature  and  pressure  ratios 
are  decreased  20  percent,  the  mass  loss  ratio  Increases  3  percent,  and  the 
density  remains  virtually  the  same.  This  Is  the  largest  effect  seen  thus  far 
In  the  sensitivity  studies,  and  must  be  kept  In  mind  when  the  comparisons  of 
experimental  data  with  predictions  are  made. 
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Fig.  5.  Sensitivity  to  Chemical  Species  Relaxation  Rates 
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Most  of  the  carbon  vapor  Is  made  up  of  molecules.®  Therefore, 
the  relaxation  of  internal  energy  states  across  the  Knudsen  layer  is  likely  to 
be  somewhere  between  frozen  and  equilibrium.  That  is,  the  internal  energy  of 
the  C3  molecule  at  moderate  temperature  {<4500*10  is  distributed  mainly 
between  rotational  and  vibrational  modes.  Since  the  rotational  energy  states 
equilibrate  (relax)  rapidly  within  a  few  molecular  collisions,  they  are  likely 
to  be  in  equilibrium.  In  contrast  to  this,  the  vibrational  energy  states 
equilibrate  (relax)  slowly  requiring  thousands  of  molecular  collisions.  Thus, 
they  are  likely  to  remain  nearly  frozen. 

d.  Surface  Adsorbed  Layer  Kinetics 


The  baseline  calculations  assume  s  a  ,  i.e.,  all  of  the  species 
vaporization  coefficients  are  equal  to  one  another.  If,- instead  of  this,  the 
Individual  species  remain  in  the  adsorbed  layer  on  the  surface  long  enough  for 
chemical  equilibrium  to  be  established  within  the  layer,  then  the  ot^  will 
not  be  equal  but  rather  related  to  through  Eq.  (8b).  The  sensitivity  of 
calculated  results  to  the  carbon  species  vaporization  coefficients  Is  shown  in 
Fig.  7. 


Results  are  shown  for  Cases  A,  8,  and  C  with  species  vaporization 
coefficients  as  given  below: 


Case 

al 

a3 

«1 

A  (baseline) 

0.39 

0.39 

«i  = 

B 

0.73 

0.39 

a  -  (c^)1 

C 

1.0 

0.39 

a1,1*3  *  1,0 

Even  with  the  extreme  (unphysical)  assumption  C,  the  calculated  results  differ 
very  little  from  the  baseline  results  (A).  This  Indicates  the  dominant 
Importance  of  cty 
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Sensitivity  to  Adsorbed  Kinetics 
(For  case  definitions,  see  text) 


B.  MOLECULAR  DYNAMICS  SIMULATIONS 


The  ad  hoc  manner  In  which  the  mass,  momentum,  and  energy  conserva¬ 
tion  equations  discussed  In  Subsection  I I -A  must  be  closed  fortunately  does 
not  lead  to  great  uncertainty  In  the  primary  calculated  variables  of  Interest. 
However,  the  conservation  equation  approach  Is,  effectively,  only  able  to  con¬ 
sider  limiting  cases  of  frozen  or  equilibrium  chemical  reaction  rates  for 
adsorbed  layer.  Internal  energy  states,  and  carbon  species  relaxation  proces¬ 
ses.  The  sensitivity  studies  indicate  that  model  predictions  for  use  In  the 
determination  of  carbon  total  vapor  pressure  and  vaporization  coefficients  are 
also  not  critically  dependent  upon  the  rate  of  these  relaxation  processes. 
Therefore,  for  the  principal  objectives  of  this  work,  a  model  based  upon  the 
conservation  equation  approach  appears  to  be^ adequate. 

If  a  Knudsen  layer  model  capable  of  Including  the  effects  of  the  above 

finite  rate  relaxation  processes  were  available,  it  should  be  possible  to 

determine  these  rates  experimentally  by  comparing  data  with  model  predictions. 

In  this  way,  a  Knudsen  layer  can  be  considered  a  tool  for  chemical  kinetics 

studies  In  a  manner  analogous  to  the  well-known  use  of  shock  waves  for  such 
1 5 

studies.  The  differences  In  these  approaches  are  appealing  from  an  instru¬ 
mentation  point  of  view.  Whereas  with  a  shock  wave  the  kinetic  processes  are 
relaxing  to  a  higher  energy  state  (temperature),  In  a  Knudsen  l^yer,  they  are 
relaxing  to  a  lower  energy  state  tnus  aiding  the  early  time  (Initial) 
detection  accuracy. 


15Resler,  E.  L..,  "The  Shock  Tube  and  Chemical  Kinetics,"  Fluid  Dynamics  and 
Applied  Mechanics,  Ed.  Diaz  and  Pal,  1961,  pp.  125-145. 


A  multi-species  chemically  reacting  Knudsen  layer  model  is  presently 
being  developed  through  a  collaborative  effort  with  the  University  of  Texas  at 
Austin.  The  model  1$  based  upon  reactive  molecular  dynamics  simulation  algo¬ 
rithms^®*^  which  have  previously  been  applied  to  complex  problems  of  non- 
equilibrium  chemistry.  Such  a  model  seeks  to  retain  all  essential  features  of 
physical  and  chemical  reality.  A  report  on  this  work  will  be  written  separately. 


16Turner,  J.  S.,  "Discrete  Simulation  Methods  for  Chemical  Kinetics,"  Journal 
_of  Physical  Chemistry,  Vol.  81,  1977,  pp.  2379-2408. 

'7Turner,  J.  s.,  "From  Microphysics  to  Macrochemistry  via  Discrete  Simula¬ 
tions,"  Computer  Modeling  of  Matter,  ACS  Symposium  Series  Mo.  86,  American 
Chemical  Society,  Washington,  b.d.,  Ed.  P.  Lykos,  pp.  231-264. 
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III.  DATA  INTERPRETATION 


A.  SUMMARY  OF  APPROACH,  CRITICAL  ASSUMPTIONS 


1.  Relationship  of  Directly  Measured  Quantities  to  Desired  Thermochemical 

Properties 

The  way  in  which  the  continuous  wave  (CW)  laser  and  the  pulsed  laser 
carbon  ablation  experiments  are  related  to  overall  program  objectives  Is  shown 
schematically  in  Fig.  8.  Whereas  both  sets  of  experimental  data  supply  vapori¬ 
zation  coefficient  information,  carbon  vapor  pressure  values  can  only  be  de¬ 
rived  from  the  CW  laser  data.  The  reason  for  this  limitation  of  information 
available  from  the  data  Is  discussed  in  Appendix  C.  Only  the  pulsed  laser  data 
reach  temperatures  sufficiently  high  to  potentially  suppTy  carbon  melting  data. 

a.  Vaporization  Coefficients 

The  present  method  of  obtaining  vaporization  coefficients  from  the 
laser  ablation  data  is  closely  related  to  the  classical  Langmuir  Method.  The 
Langmuir  vapor  pressure  pL  is  the  product  of  the  vaporization  coefficient 
and  the  equilibrium  vapor  pressure  p$  of  the  saturated  vapor  at  temperature 
Ts.  Classically,  these  quantities  are  related  to  the  vaporization  rate  m 
and  the  molecular  weight  M  by  the  Langmuir  relation 

/27TRT 

PL  »  aps  =•  my — (10) 

The  Langmuir  method  Is  carried  out  by  vaporizing  material  into  a 
vacuum.  The  analysis  discussed  In  Subsection  II-A  Indicates  that  when  vapori¬ 
zation  takes  place  Into  a  finite  back  pressure  environment,  but  the  equilibrium 
v.apo,r  pressure  p$  Is  sufficiently  greater  than  the  ambient  pressure  p#. 
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the  Langmuir  pressure  and  the  vaporization  (mass  loss)  rate  m  are  related  by 


I 2  jtrt 

PL  *  aps  *  F^o.s)  my— — 
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where  Fj(a,  u)  is  a  nonlinear  function  of  vaporization  coefficient  a  and 
also  depends  on  the  ratio  of  specific  heats  8. 

Thus,  If  Eq.  (11)  Is  used  Instead  of  the  Langmuir  relation,  Eq.  (10), 
measured  mass  loss  rate  m  data  can  be  converted  to  Langmuir  pressures.  Then 
for  a  given  vapor  pressure  relationship  which  supplies  ps  as  a  function  of 
temperature  T  ,  a(T$)  can  be  determined.  Because  the  function  F^a,8) 
is  nonlinear,  the  conversion  of  m  to  a  is  best  accomplished  by  carrying  out 
the  Iteration  graphically.  This  is  Illustrated  in  Subsection  II I -B.  For  the 
NASA  pulsed  laser  experiments,  the  Langmuir  pressure  pL  is  similarly  related 
to  the  measured  free-jet  stagnation  pressure  pQ,  i.e. 

pL  =  F2(<*,#)p0  (12) 

It  should  be  emphasized  that  the  values  of  vaporization  coefficient 
as  a  function  of  temperature  determined  in  this  manner  depend  on  the  equili¬ 
brium  vapor  pressure,  i.e.,  a  =  pL/p$.  To  uniquely  determine  the  vapori¬ 
zation  coefficient  from  pL,  the  correct  equilibria  vapor  pressure  relation¬ 
ship  must  be  known  as  discussed  below. 

b.  Vapor  Pressure 

Once  the  Langmuir  vapor  pressure  as  a  function  of  Ts  Is  known, 
vaporization  coefficients  aA(T$)  and  <*g(Ts)  corresponding  to  literature  vapor 
pressure  relationships  A  and  B,  respectively,  can  be  determined.  The  analy¬ 
tical  model  discussed  In  Subsection  II-A  can  then  be  used  to  predict  m  as  a 


function  of  T$  for  constant  p#.  Calculations  of  this  type  for  two 
representative  carbon  vapor  pressure  specifications,  A  and  B  (not  to  be 
confused  with  Cases  A,  B,  and  C  of  Section  II),  are  shown  In  Fig.  9  for 
ambient  pressures  pa  of  1  atm  and  0.1  atm.  For  a  given  pa,  the  predicted 
results  are  asymptotic  to  Eq.  (11)  at  the  upper  end  of  the  temperature  range, 
l.e.,  trfiere  p. »p.  .  Here,  both  thermochemical  data  sets  A  and  B,  with 
vaporization  coefficients  aA(Ts)  end  <*B(TS),  respectively,  predict 
the  same  values  of  total  mass  loss  rate  A. 

For  given  pa,  as  the  temperature  T  Is  decreased  the  predicted  M 
values  for  vapor  pressure  relationships  A  and  B  differ.  This  difference 
becomes  greater  as  T$  Is  lowered  still  more.  At  the  lower  range  of  tempera¬ 
tures  for  given  p#,  H  varies  rapidly  with  T$  and  the  predicted  results 
become  asymptotic  to  the  linear  Knudsen-langmul r  equation 

*  ■  ay[z&rs  "3> 

where  Asp  Is  the  area  of  the  laser  Irradiated  spot. 

From  Eq.  (13),  It  can  easily  be  seen  that  the  rapid  Increase  of  M  with 
T$  occurs  whenever  PS(TS)  Is  only  slightly  greater  than  pa.  Since  vapor 
pressure  relationships  A  and  B  predict  p.  »  p,  at  different  values  of  T  , 

t  5  d  5 

the  predicted  differences  In  M  versus  1/T$  behavior  seen  In  Fig.  9  mqy  be 
directly  attributed  to  the  different  values  of  pg  at  a  given  T  . 

Given  measurements  of  M  and  T.  at  a  specified  value  of  pa,  com- 
pari  son  of  such  experimental  data  with  predicted  results  such  as  those  shown  In 
Fig.  9,  allows  the  proper  equilibrium  vapor  pressure  relationship  (A  or  B  or 
some  other)  to  be  determined. 
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c.  Melt  Temperature 


As  originally  envisioned  in  this  work,  the  occurrence  of  carbon  melt¬ 
ing  was  to  be  related  to  observed  changes  in  M  versus  1/T$  behavior  of  the 
high  temperature  pulsed  laser  data.  Since  the  heat  of  vaporization  (liquid- 
vapor)  Is  always  less  than  the  heat  of  sublimation  (sol id- vapor),  an  observed 
change  of  slope  In  the  data  could  possibly  be  related  to  the  onset  of  melting. 

Subsequently,  a  more  direct  approach  Involving  monitoring  the  surface 

1*1 

reflectivity  as  a  function  of  time  during  the  laser  pulse  was  recommended 

and  Implemented.  For  the  laser  annealing  of  semiconductor  materials,  this 

technique  has  been  used  very  successfully  to  detect  melting  by  observing  the 

18 

large  Increase  In  reflectivity  when  the  semiconductor  surface  melts.  An 
auxiliary  laser  is  used  to  monitor  the  reflectivity. 

2.  Critical  Assumptions 

The  purpose  of  the  following  discussion  Is  to  examine  some  of  the  key 
assumptions  related  to  the  validity  of  the  data  interpretations  presented  in 
Subsection  III-B. 

a.  Continuous  Wave  Laser  Ablation  Data 

For  this  data  set,  the  directly  measured  quantities  are  total  mass 
loss,  run  time,  and  surface  temperature.  The  quantities  compared  to  analytical 
model  predictions  are  the  mass  loss  rate  M  and  the  saturated  vapor  temperature 
T$  which  Is  assumed  to  be  equal  to  the  measured  surface  temperature. 


l®Au$ton,  D.  H.,  et  al.,  "Time  Resolved  Reflectivity  of  Ion-Implanted  Silicon 
During  Laser  Annealing,"  Applied  Physics  Letters,  Vol.  33,  1978,  pp.  437-440. 
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To  obtain  the  mass  loss  rate,  two  key  assumptions  are  made.  First,  It 
Is  assumed  that  all  mass  loss  takes  place  by  sublimation  and,  secondly,  that 
the  area  over  which  the  sublimation  takes  place  Is  the  same  as  the  Irradiating 
laser  spot  size,  even  though  the  sample  Is  rotating.  It  Is  difficult  to  con¬ 
clusively  substantiate  these  assumptions.  The  latter  assumption  Is  at  least 
qualitatively  confirmed  by  the  visually  observed  carbon  vapor  plume,  l.e.,  its 
flow  direction  Is  nearly  parallel  to  the  axis  of  the  Incoming  laser  beam.  The 
differential  time  used  to  convert  the  total  mass  loss  to  a  mass  loss  rate  was 
obtained  from  the  strip  chart  records  of  surface  temperature.  The  time  Incre¬ 
ment  used  in  each  case  was  that  for  which  the  surface  temperature  remained 
within  a  small  fraction,  approximately  5  percent  of  its  highest  value  T$. 

The  most  critical  assumption  related  to  the  measured  surface  tempera¬ 
tures  is  the  assumed  emisslvity  value  of  unity.  For  the  pyrolytic  graphite 

samples  which  were  irradiated  and  viewed  by  the  pyrometer  in  a  direction 
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parallel  to  the  a-b  plane,  this  assumption  Is  most  defensible.  *  A  change 
of  emlssivity  from  1.0  to  0.9  shifts  the  measured  surface  temperatures  up 
~  100°K  which  would  significantly  affect  the  vapor  pressure  conclusions. 

b.  Pulsed  Laser  Ablation  Data 


For  this  data  set,  the  stagnation  pressure  pQ  of  the  vaporized 

carbon  plume  Is  assumed  to  be  derivable  from  the  measured  Mach  disk  geometry 

21 

using  a  simple  algebraic  relationship.  This  relationship  Is  substantiated 


19Pyrolyt1c  Graphite  Handbook,  Metallurgical  Production  Dept.,  General  Electric 
co.,  Schenectady,  HT,  1964,  as  quoted  In  Gokcen,  N.  A.,  et  al.,  "Determination 
of  Graphite/Ll quId/Yapor  Triple  Point  by  Laser  Heating,"  High  Temperature 
^clence,  Vol.  8,  June  1976,  pp.  81-97. 

z°foulouk1an,  Y.  $.,  Ed.,  Thermophysical  Properties  of  Matter,  Vol.  8,  IFI/ 
Plenum  Press,  NY,  1970. 

z'Cov1ngton,  M.  A.,  G.  N.  Liu,  and  K.  A.  Lincoln,  "Free-jet  Expansions  from 
Laser-vaporized  Planar  Surfaces,"  AIAA  Journal,  Vol.  15,  August  1977,  pp.  1174- 
1179.  - 
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for  sonic  orifices  for  a  wide  range  of  gases,  temperatures,  and  pressure 
22 

ratios.  It  Is  further  assumed  that  the  choked  carbon  vapor  flow  vaporized 
from  the  solid  carbon  surface  subsequently  behaves  as  If  emanating  from  an 
orifice.  The  pyrometrlcally  measured  surface  temperature,  as  In  the  CW  laser 
data  case.  Is  sensitive  to  the  assumed  emlsslvlty. 

The  relatively  steep  slope  of  stagnation  pressure  pQ  versus  1/T 
observed  In  the  pyrolytic  graphite  vaporization  data  necessitated  considerable 
reexamination  of  the  methods  for  determining  both  p0  and  T$.  This  slope 
was  expected  to  have  a  numerical  value  near  that  of  the  normalized  heat  of 
formation  of  carbon,  l.e.,  4H*/R  «  93,000*K.  The  slope  seen  in  the  data  is 
about  126,000*K.  Possible  explanations  for  this  behavior  which  have  been 
considered  are: 

o  Emlsslvlty  a  function  of  temperature 
o  Steady-state  energy  balance  not  achieved 
o  Significant  radial  conduction  of  energy 
o  Unsteady  surface  morphology 

o  Condensation  of  carbon  vapor  In  the  free- jet  plume 

o  Energy  release  effects  In  the  fnee-jet  plume 

o  Strongly  temperature-dependent  vaporization  coefficients 

Our  conclusions  thus  far  regarding  these  effects  are  discussed  consecutively 
below. 


If  the  emlsslvlty  changes  with  temperature,  it  Is  likely  to  Increase 
as  temperature  Increases.  It  can  easily  be  shown  that  accounting  for  such  a 
temperature  dependence  In  the  data  reduction  would  increase  the  slope  of  the 
measured  data,  making  It  even  more  difficult  to  explain. 


22Ashkenas,  H.  and  F.  S.  Sherman,  Rarefied  Gas  Dynamics,  Ed.  J.  H.  de  Leeuw, 
Academic  Press,  NY,  1966,  Vol.  2,  p.  64. 
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The  effect  of  not  having  achieved  a  steady-state  balance  of  energy  at 
the  ablating  surface  has  been  looked  at  by  performing  one-dimensional  heat  con¬ 
duction  calculations.  Even  at  the  lowest  laser  Input  energy  levels,  where  the 
time  to  steady  state  would  be  longest,  a  reasonably  quasi-steady  Mach  disk  flow 
field  Is  established.  Thus,  non-steady-state  energy  balance  effects  do  not 
appear  to  be  Influencing  the  measured  data. 

For  the  pulsed  laser  data,  the  pyrolytic  graphite  samples  were  Irradi¬ 
ated  parallel  to  the  c-axis.  Thus,  the  high  thermal  conductivity  direction  was 
radially  outward  from  the  circular  laser  spot.  To  determine  if  radial  conduc¬ 
tion  of  energy  was  significantly  affecting  the  measured  data,  additional  pyro¬ 
lytic  graphite  data  were  obtained  using  a  new  sample  configuration  to  minimize 
radial  conduction.  Cylindrical  samples  were  Irradiated  end-on,  with  the 
cylinder  diameter  approximately  equal  to  the  laser  spot  size.  The  new  data 
showed  no  significant  differences  from  the  older  data  for  which  the  laser  spot 
was  focussed  on  a  flat  planar  area  much  larger  than  the  spot  size.  Thus,  It 
was  concluded  that  radial  conduction  of  energy  was  not  affecting  the  data. 

As  ablation  (surface  mass  removal)  Is  Initiated,  the  surface  morpho¬ 
logy  changes  as  etch  pits  are  formed  and  the  number  of  sites  at  which  atoms 
and  molecules  are  released  from  the  solid  increases.  When  the  surface  has 
been  completely  covered  with  etch  pits,  the  nimber  of  sites  releasing  atoms/ 
molecules  remains  fixed,  i.e. ,  the  surface  morphology  is  at  steady  state.  At 
the  lower  levels  of  laser  Irradiation,  considerably  less  mass  Is  removed  from 
the  sample  than  at  the  highest  level  of  energy  deposition.  If  the  amount  of 
mass  removed  Is  not  sufficient  to  have  established  a  steady-state  surface, 
then  Increasing  numbers  of  surface  vaporization  sites  with  Increasing  energy 
deposition  could  be  the  cause  of  the  unexpectedly  high  slope  In  the  data. 

Based  upon  estimates  of  the  amount  of  mass  removal  required  to  establish  a 

13  23 

steady-state  surface  morphology,  ’  It  was  concluded  that  a  steady  surface 
Is  likely  to  be  established  even  at  the  lowest  levels  of  energy  deposition. 
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The  analytical  calculations  Indicate  that  If  the  chemical  reaction 
rates,  l.e. ,  the  carbon  species  equilibration  rates  as  well  as  the  Internal 
energy-translational  energy  relaxation  rates,  are  frozen,  then  the  expanding 
carbon  vapor  leaving  the  nonequilibrium  region  close  to  the  surface  (the 
Knudsen  layer)  Is  highly  supersaturated.  If  this  vapor  subsequently  condenses 
to  solid  carbon  In  the  free-jet  expansion,  measured  stagnation  pressures  would 
be  too  low.  To  "explain11  the  data  slope,  a  physical  explanation  supporting 
significant  condensation  at  the  lower  end  of  the  temperature-stagnation  pres¬ 
sure  range  and  little  or  no  condensation  at  the  upper  end  would  be  required. 

It  Is  planned  In  future  experiments  to  look  for  the  presence  of  condensed 
carbon  by  measuring  backscatterlng  of  He-Ne  laser  light  from  the  free-jet  plume 
region. 


There  are  two  potentially  significant  sources  of  energy  release  In  the 
expanding  carbon  vapor  free-jet.  A  frozen  vapor  at  the  edge  of  the  Knudsen 
layer  would  begin  to  react  and  equilibrate  carbon  species  In  the  free-jet 
region  thereby  releasing  about  200  cal/mol.  Similarly,  nonequilibrium  between 
Internal  and  translational  energy  degrees  of  freedom  at  the  Knudsen  layer  edge 
would  also,  upon  equilibration,  release  about  the  same  amount  of  energy.  Argu¬ 
ing  that  this  energy  Is  transferred  Into  translational  degrees  of  freedom  and 
thus  Increases  the  pressure  (and  the  stagnation  pressure),  the  measured  stagna¬ 
tion  pressures  would  be  too  high  whenever  significant  colli slonal  relaxation 
takes  place.  Since  colllslonal  relaxation  increases  with  pressure  and  the 
measured  stagnation  pressures  are  "too  high"  at  larger  values  of  pressure,  this 
effect  could  explain  the  data.  However,  literature  data  for  Mach  disk  geome¬ 
tries,  including  energy  release  effects  (relaxation  of  internal  and  transla¬ 
tional  degrees  of  freedom),  shows  a  very  weak  dependence  on  the  ratio  of 
24 

specific  heats  8.  This  Indicates  that  energy  release  has  little  effect  on 
the  Mach  disk  location  (geometry). 


^Rosenblatt,  G.  M.,  "The  Role  of  Defects  In  Vaporization:  Arsenic  and  Anti¬ 
mony,"  Surface  Defect  Properties  of  Solids,  Yol.  5,  1976,  pp.  36-64. 

2*Wen,  Y.  S.,  "derivation  of  the  Free-jet  Mach-disk  Location  Using  the  Entropy- 
balance  Principle,"  Physics  of  Fluids,  Vol.  18,  No.  11,  pp.  1421-1425. 
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If  the  vaporization  coefficient  changes  rapidly  with  temperature ,  the 

temperature  dependence,  which  Is  expressed  through  an  activation  energy  4H*, 
25 

Is  represented  by 


in  a  =  +  4S*  -  AS #  (14) 

where  AS *  Is  an  activation  entropy.  Combining  this  with  the  vapor  pressure 
expression  which  is  given  in  terms  of  the  standard  heat  of  formation  AH *  and 
the  standard  entropy  of  formation  by 


in  Ps  =  (15) 

the  Langmuir  vapor  pressure  pL  =  otp$  is  given  by 


in  PL  =  +  4S*  (16) 

From  Eq.  (14),  a  rapid  variation  of  a  with  T  requires  4H*  >>  4H*.  Combining 
Eqs.  (12)  and  (16) 


in  P0  =  -ffr —  +  ^S*-in  F2(«,») 


(17) 


25 


Rosenblatt,  G.  M.,  "Evaporation  from  Solids,"  Treatise  on  Solid  State 
Chemistry,  Vol.  6a,  Ed.  N.  B.  Hannqy,  Plenum  Press,  kY,  I&/6,  p.  199. 
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The  dominant  importance  of  the  numerical  value  of  the  carbon  species 

vaporization  coefficient  a ^  on  the  calculated  results  was  shown  by  the 

parametric  sensitivity  study  discussed  In  relation  to  Fig.  7  In  Subsection 

II-A.  Because  of  this  and  because  the  assumption  of  equal  values  of  vaporlza- 

1 3 

tlon  coefficients  for  all  species  Is  reasonable  above  3000 *K,  In  all  the 

results  presented  here  It  has  been  assumed  that  s  a. 

For  given  values  a  and  8  ,  the  calculated  results  shown  In  Fig.  3 
and  represented  functionally  by  Eq.  (11),  can  be  expressed  as 


m  *  c^y^-Pj  (18) 

where  C  is  a  constant.  A  semi-log  plot  of  m  versus  1/T$  [see  Eq.  (15)]  thus 
will  have  nearly  the  same  slope  as  ps,  l.e. ,  the  slope  will  be  -4H*/R  where 
4H*  Is  the  reference  heat  of  formation.  For  different  values  of  a,  the 
"constant"  C  has  a  different  numerical  value.  Lines  showing  M  =»  mASp 
versus  1/T$  for  the  various  values  of  a  will  be  parallel  to  one  another  as 
shown  In  Fig.  10.  The  total  mass  loss  rate  M  (mg/sec)  Is  obtained  from 
m  (g/cn£-sec)  by  multiplying  by  the  laser  spot  area  A^fcrrt2)  and  dividing  by 
1000. 
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The  lines  in  Fig.  10  were  calculated  assuming  the  total  vapor  pressure 
Ps  to  be  given  by  the  carbon  species  thermochemlcal  data  tabulated  in  Refs. 

26  and  27.  Hereafter,  these  thermochemical  data  will  be  referred  to  as  the 
Livermore  data.  If,  Instead,  the  JANAF  data^  are  used,  the  predicted  M  versus 
inverse  temperature  behavior  will  be  as  shown  in  Fig.  11.  Some  of  the  reasons 
for  the  difference  in  these  two  thermochemical  data  sets  are  discussed  in  the 
final  paragraphs  of  Subsection  II1-B-2.  Note  the  different  slopes  and  inter¬ 
cepts  for  given  constant  a  in  Figs.  10  and  11. 

Comparison  of  the  predicted  results  in  Figs.  10  and  11  with  experi¬ 
mental  data  allows  a  to  be  determined  graphically  as  a  function  of  T  .  Low 
ambient  pressure  data  for  pyrolytic  graphite  samples  obtained  in  the  continuous 
wave  laser  experiments  (See  Appendix  A  and  Ref.  10)  are  also  shown  in  these 
figures  along  with  a  best  fit  line  through  the  data.  Comparison  of  the  data 
lines  in  Figs.  10  and  11  with  the  predicted  results  gives  the  vaporization 
coefficient  values  and  temperature  dependence  shown  in  Fig.  12  for  the  two 
thermochemlcal  data  sets. 

Similarly,  for  given  values  of  a  and  8,  the  calculated  results  shown 
In  Fig.  13  and  represented  functionally  by  Eq.  (12),  can  be  expressed  as 

P0  *  C2ps  (19) 


26Lee,  E.  L.  and  R.  H.  Sanborn,  "Extended  and  Improved  Thermal  Functions  for 
the  Gaseous  Carbon  Species  C1-C7  from  298  to  10,000'K,"  Hiqh  Temperature 

Science,  Vol.  5,  1973,  pp.  438-453.  — * - - 

2 /L elder,"  H.  R.,  0.  H.  Krlkorlan,  and  D.  A.  Young,  "Thermodynamic  Properties  of 
Carbon  Up  to  the  Critical  Point,"  Carbon,  Yol.  11,  1973,  pp.  555*563. 
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Fig. 


Comparison  of  CW  Laser  Data  and  Predicted 
Results  (JANAF) 
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Fig.  12.  Vaporization  Coefficients  Derived  from  CM  Laser  Data 


Figures  13  and  14  give  the  predicted  stagnation  pressure  pQ  as  a  function  of 
Inverse  temperature  for  the  same  thermochemlcal  data  and  vaporization  coeffi¬ 
cient  values  as  In  Figs.  10  and  11.  The  pulsed  laser  stagnation  pressure  data 
are  presented  In  Fig.  15.  The  data  line  from  Fig.  15  Is  also  shown  In  Figs.  13 
and  14.  (It  Is  severely  truncated  In  Fig.  14  because  a  values  >1.0  are  not 
physically  meaningful.)  Vaporization  coefficients  obtained  grahlcally  from 
Figs.  13  and  14  are  seen  In  Fig.  16. 

The  results  of  Figs.  12  and  16  are  combined,  and  the  experimentally 
determined  values  of  a  from  3400 *K  to  4600 *K  are  presented  in  Figs.  17  and  18 
for  the  Livermore  and  JANAF  thermochemlcal  data,  respectively.  In  Fig.  17,  the 
high  sensitivity  of  the  results  to  the  value  of  surface  emlssivity,  used  in 
reducing  the  pyrometer  data  to  obtain  surface  temperature.  Is  also  shown.  An 
emlssivity  of  0.97  for  the  pulsed  laser  data  Is  particularly  attractive  because 
the  numerical  values  of  the  high  temperature  data  then  match  the  lower  tempera¬ 
ture  data  at  4000 *K  and  approach  unity  at  the  Livermore  theoretical  triple 
point  (melt)  temperature  (4765*K).  In  marked  contrast,  the  vaporization 
coefficients  based  on  the  JANAF  thermochemical  data  (Fig.  18)  are  physically 
meaningful,  l.e.,  less  than  unity,  above  4200#K  only  If  much  lower  values  of 
emlssivity  are  appropriate  for  the  pulsed  laser  data. 

2.  Vapor  Pressure 

The  vaporization  coefficient  discussion  just  completed  gives  values 
of  the  vaporization  coefficient  as  a  function  of  temperature  derived  from  two 
vapor  pressure  relationships,  i.e.,  Livermore  and  JANAF.  The  total  vapor 
pressure  for  carbon  predicted  by  these  two  sets  of  recommended  thermochemlcal 
data  is  shown  in  Fig.  19.  The  Livermore-predicted  total  vapor  pressure  Is 
higher  by  a  factor  of  about  2.5  at  4000 *K  and  by  a  factor  of  about  3.9  at 
4765*K . 
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Fig.  13.  Comparison  of  Pulsed  Laser  Data  and  Predicted 
Results  (Livermore) 
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Fig.  17.  Vaporization  Coefficients  for  Livermore  Thermochemical  Data 
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Fig.  18.  Vaporization  Coefficients  for  JANAF  Thermochemical  Data 


The  intended  method  of  determining  which  thermochemical  data  are  more 
nearly  correct  has  been  reviewed  in  Subsection  III-A  and  is  discussed  in  detail 
in  Ref.  10.  It  is  based  upon  predicted  differences  in  behavior  between  experi¬ 
mental  data  obtained  when  the  ambient  pressure  is  low  relative  to  the  vapor 
pressure  and  other  experimental  data  obtained  for  higher  ambient  pressures  (the 
order  of  the  vapor  pressure).  This  is  illustrated  in  Fig.  20. 

The  data  points  used  to  determine  the  temperature  dependence  of  the 
vaporization  coefficients  and  the  consequent  “agreement"  of  JANAF-  and 
Livermore-predicted  mass  loss  rates  with  the  data  are  seen  in  Fig.  20a.  The 
same  data  points  and  predicted  mass  loss  rates  for  ambient  pressure  p.  = 

a 

1  Torr  are  shown  in  Figs.  20b  and  20c.  Also  shown  in  these  latter  figures  are 
predicted  mass  loss  rates  for  p.  ranging  from  0.02  to  1.0  atmospheres  for  the 

a 

JANAF  and  Livermore  thermochemical  data  with  corresponding  vaporization  coeffi¬ 
cients.  In  Figs.  21  through  23,  these  predicted  results  are  compared  with  the 
experimental  data. 

The  data  are  very  sparse,  especially  in  Fig.  21.  For  these  data,  the 
ambient  pressure  is  not  high  enough  relative  to  the  vapor  pressure  to  cause  a 
measurable  difference  in  JANAF-  and  Livermore-predicted  results  at  the  measured 
surface  temperature  T$.  For  the  higher  ambient  pressure  data  in  Figs.  22  and 
23,  the  observed  dependence  of  mass  loss  rate  on  ambient  pressure  and  surface 
temperature  is  reasonably  in  agreement  with  analytical  predictions  based  upon 
the  Livermore  thermochemical  data.  In  contrast  to  this,  much  less  satisfactory 
agreement  of  the  experimental  data  with  predictions  is  obtained  when  JANAF 
vapor  pressures  are  used.  Based  upon  this  and  the  discussion  below,  it  is 
concluded  that  the  Livermore  data  should  be  used  to  calculate  carbon  vapor 
pressure  as  a  function  of  temperature. 

The  basis  for  the  Livermore  thermochemical  data  appears  to  be  very 
well  founded.  The  reference  heat  of  formation  used  for  the  critical  carbon 
species  Cj  molecule  is  within  0.5  kcal/mol  of  the  average  experimental  value 
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Fig.  20.  Variation  of  Mass  Loss  Rate  with  Surface  Temperature  and  Ambient  Pressure 
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determined  by  Wachi  and  Gilmartln.  This  Is  significant  because  these 

Investigators,  for  the  first  time,  obtained  very  good  agreement  of  their  third- 

law  and  second-law  derived  heat  of  formation.  The  free  energy  functions  for 

26 

the  Livermore  data  were  calculated  by  Lee  and  Sanborn  Including  contribu¬ 
tions  of  electronic  degrees  of  freedom.  Without  the  very  high  temperature 
contributions  due  to  the  electronic  degrees  of  freedom,  their  thermal  functions 
are  essentially  in  agreement  with  the  established  values  of  Strauss  and 

Thiele,  and  their  entropy  values  are  above  the  "lower  limits"  determined 

30 

by  Hansen  and  Pearson.  In  contrast  to  this,  the  JANAF-recomnended  entro¬ 
pies  for  are  below  the  Hansen  and  Pearson  lower  limits,  and  the  JANAF 
reference  heat  of  formation  for  C3  Is  5.7  kcal/mol  less  than  the  average 
Wachi  and  Gllmartin  value. 

A  final  appealing  feature  of  the  Livermore  data,  as  discussed  by 
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Leider,  Krikorian,  and  Young,  is  that  the  total  vapor  pressure  reaches  the 
triple  point  pressure  (wlOO  atm)  at  4765*K.  This  is  several  hundred  degrees 
above  the  literature  values  for  triple  point  (melt)  temperature,  but  consistent 
with  the  melt  temperature  conclusion  discussed  In  the  following  section.  A 
significant  fraction  of  the  Livermore-predicted  total  carbon  vapor  pressure  at 
4765*K  is  due  to  carbon  species  above  Cg,  especially  C7.  The  contributions 
of  these  species  to  the  total  vapor  pressure  are  supported  by  the  larger  slope 
of  the  high  temperature  Langmuir  vapor  pressure  data  in  Fig.  15.  The  JANAF 
tables  presently  do  not  give  data  for  carbon  species  above  Cg.  Because  of 
this  and  the  low  JANAF  C3  vapor  pressure,  the  JANAF- predicted  triple  point 
temperature  is  «5125"K. 


2®Wach1,  F.  N.  and  D.  E.  Gllmartln,  "Heat  of  Formation  and  Entropy  of  C3  Mole- 
cule,4  High  Temperature  Science,  Vol.  4,  No.  5,  pp.  423-431. 

"Strauss,  H.  L.  and  E.  TMele,  "Thermodynamics  of  C3.  II.  General  Methods 
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Melt  Temperature 

Carbon  Phase  Change  Detection  Using  Surface  Reflectance 

This  investigation  of  the  possible  phase  change  of  carbon  used  a  new 
technique  evidently  not  employed  previously  in  carbon  studies  to  detect  the 
presence  of  a  melt  layer  on  a  solid  substrate.  The  application  of  this  method 
was  inspired  by  recent  success  in  the  use  of  time-resolved  spectral  reflecti¬ 
vity  of  doped  silicon  to  monitor  surface  phase  change  during  laser-heated 

lO 

annealing.  The  reflective  behavior  of  germanium  (Ge)  and  silicon  (Si)  heated 
to  melting  are  included  along  with  carbon  results  to  provide  comparative  data 
when  a  liquid  layer  was  known  to  be  present. 

The  reflection  and  refraction  of  an  electromagnetic  wave  at  the  optic¬ 
ally  smooth  surfaces  of  liquid  and  solid  materials  can  be  accurately  described 
by  the  Fresnel  equations  in  terms  of  the  material's  optical  indices  and  the 
incidence  angle  of  the  radiation.  However,  once  a  surface  is  heated  to  tem¬ 
peratures  where  melting  or  significant  vaporization  occurs,  it  becomes  rough¬ 
ened  and  is  no  longer  optically  smooth.  The  induced  discontinuities  cause 
scattering  of  incident  radiation  independent  of  the  material's  intrinsic 
optical  properties  and,  if  the  size  of  the  irregularities  is  of  the  order  of 
the  radiation  wavelength  or  larger,  the  scattering  dominates  the  reflective 
behavior.  In  this  case,  the  directional  distribution  of  reflected  radiation 
can  be  best  described  by  geometric  optics  that  take  into  account  the  statisti¬ 
cal  nature  of  surface  features. 
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In  these  experiments.  It  was  found  to  be  advantageous  to  roughen  the 
Ge  and  SI  surfaces  prior  to  exposure  to  Insure  uniform  melting  over  the  heated 
area.  Also,  graphite  surfaces  were  roughened  by  vaporization  regardless  of 
their  Initial  smoothness.  This  makes  the  strict  application  of  the  Fresnel 
relations  for  predicting  expected  changes  between  solid  and  liquid  phases  in¬ 
appropriate  for  the  data  included  herein.  Since  the  purpose  of  this  investi¬ 
gation  was  to  measure  relative,  transient  changes  in  surface  reflectance  as  a 
means  to  detect  melting,  comparisons  between  absolute  reflectivity  values  pre¬ 
dicted  from  electromagnetic  theory  for  smooth  surfaces  and  experimentally 
measured  bidirectional  or  spatially  integrated  reflectivities  are  not  Included. 

However,  for  the  directional  reflectance  method  and  the  materials 
used,  the  reflected  radiation  flux  is  characterized  by  a  specular  component 
and  a  diffuse  component  distributed  about  the  reflection  angle  of  the  specular 
component.  The  specular  component  Increases  In  relative  intensity  with  In¬ 
creasing  Incidence  angle  (l.e.,  larger  grazing  angles),  and  reflection  of  a 
collimated  beam  approaches  that  for  pure  specular  reflection.  This  makes  the 
Fresnel  equations  useful  for  comparing  beam  polarization  and  Incidence  angle 
effects  as  well  as  the  relative  reflectances  expected  during  heating. 

For  a  solid  or  liquid  material  having  optical  properties  that  are 

constant  In  all  directions,  the  Fresnel  equations  for  an  Isotropic,  absorbing 

medium  relate  the  fraction  of  Incident  radiation  that  Is  reflected  from  an 

optically  smooth  surface  to  the  refractive  Index,  the  absorption  Index,  and 
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the  Incidence  angle  of  the  radiation: 

b  .  Ccosl-  a)2  *  fi2 

1  (cos*  +  00*  +  (20) 

2  2 
r  .  (bcosfl  -  a)  »  (dcosg  -  a) 

"  (#cos0  +  a)2  ♦  (6cos0  +  a)‘ 


31Born,  M., 
1964. 


and  E.  Wolf,  Principles  of  Optics,  2nd  Edition,  Permagon  Press,  MY, 
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where 

Rx  and  R)(  *  the  reflectivity  of  radiation  polarized  perpendicular 

and  parallel,  respectively,  to  the  plane  of  incidence. 

|  2“*  |  .  [(n2-k2-  s1n20)2  ♦  4n2  k2]''*  ♦  („2  -  k2  -  s1n2«) 

8  =  n2  -  k2  ,  <J  *  2nk 
n  «  refractive  Index 
k  =  absorption  index 
e  =  incidence  angle 

These  relations  are  appropriate  for  many  solid  and  liquid  substances 
including  the  crystalline  and  liquid  fie  and  Si  and,  probably,  the  liquid  phase 
of  carbon  considered  herein. 

For  strongly  absorbing  anisotropic  crystals,  the  above  relationship 
for  Isotropic  materials  does  not  properly  account  for  different  optical  proper¬ 
ties  along  the  crystal  axes.  Anisotropic  uniaxial  materials  having  only  one 
unique  direction  are  typified  by  pyrolytic  graphite  with  its  stacked  parallel 
layers  of  carbon  atoms.  The  reflection  from  optically  smooth  surfaces  of  uni¬ 
axial  crystals  having  optical  properties  along  the  c-axis  that  are  different 
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from  those  normal  to  the  axis  Is  given  by 

(cose  -  ai)2  +  b2 

R  =  - y - y - 

1  (cose  +  a1)^  +  b* 

2  2 
(c.cose-  aM  )  +  (d.cos  -  b„  ) 

R  ■  - _  ——————  -  .  y 

"  (c.cose*  a„  )  *  (d.cos  ♦  b„  ) 


(22) 


(23) 


^^Mosteller,  L.  P. ,  and  F.  Wooten,  "Optical  Properties  and  Reflectance  of  Uni¬ 
axial  Absorbing  Crystals,"  J.  Optical  Society  of  America,  Vol.  58,  19$8, 
pp.  511-518. 
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where 


Ra  and  R() 


the  reflectivity  of  radiation  polarized  perpendicular 
and  parallel,  respectively,  to  the  plane  of  Incidence. 


0}-  €2 


id  i  •  f 


(nz  '  kZ  ■  Slr,2&)2  *  4n2  kz]1/2  1  {n2  "  kZ  '  Si"2e) 


(nx"z  '  kxkz} 


(nxkz  -  n2kx) 


V 
6  • 


refractive  Index  In  j-dl recti on 
absorption  Index  In  j-dlrection 
Incidence  angle 


Note  that  the  equations  for  uniaxial  crystals  differ  from  those  for 
Isotropic  materials  only  by  Including  the  differences  In  refraction  and 
absorption  Indices  along  orthogonal  axes. 


Using  Eqs.  (20)  and  (21)  with  optical  properties  taken  from  Refs.  33 
and  34,  for  the  reflection  from  optically  smooth  surfaces  of  solid  and  liquid 
SI  In  vacuum  as  shown  in  Fig.  24,  It  is  easily  seen  that  the  use  of  a  colli¬ 
mated  beam  polarized  so  that  its  electric  vector  is  perpendicular  to  the  plane 
of  Incidence,  and  large  grazing  angles  will  maximize  the  reflected  energy. 
Likewise,  the  use  of  Eqs.  (22)  and  (23)  for  the  reflection  from  smooth  sur¬ 
faces  of  solid  pyrolytic  graphite  with  Its  anisotropic  optical  properties 
taken  from  Ref.  35  provides  a  similar  conclusion  as  shown  In  Fig.  25.  A  pre¬ 
diction  for  the  reflectivity  of  liquid  carbon  is  not  Included  since  the  opti¬ 
cal  properties  of  the  liquid  phase  are  not  known,  and  the  exact  structure  of 
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the  liquid  Is  still  a  matter  of  debate.  *  However,  there  Is  no  reason  to 
expect  that  the  reflectance  from  liquid  carbon  would  be  lower  than  that  for 
the  solid. 

b.  Experimental  Results  and  Conclusions 

A  schematic  of  the  experimental  apparatus  utilized  is  shown  in 
Fig.  26.  Material  samples  were  heated  within  a  cylindrical  vacuum  tank  having 
windows  on  both  ends  and  around  the  perimeter.  Heating  was  provided  by  a 
pulsed  Nd/glass  laser  producing  a  uniform  cross-ssct ion  beam  at  1.06  fim 
wavelength  and  operated  in  the  normal,  or  burst,  node.  The  Nd/glass  laser 
beam  irradiated  the  material  sample  normal  to  the  front  surface  over  heated 
areas  2.5  to  4  ran  in  diameter,  and  the  heating  pulse  was  approximately  tri¬ 
angular  in  shape  with  a  total  pulse  length  of  about  1.2  msec.  During  laser 
heating  cycles  the  tank  was  filled  with  argon  at  pressures  at  10  Torr  or  less. 
The  incident  beam  for  reflectivity  measurements  was  obtained  from  a  He/Ne  laser 
focused  onto  the  heated  sample  area  at  a  grazing  angle  of  60s  from  the  normal. 

The  He/Ne  laser  beam  was  polarized  perpendicular  to  the  plane  of  inci¬ 
dence  (E-vector  parallel  to  surface)  to  maximize  the  reflectance  as  previously 
noted,  and  the  beam  was  chopped  at  a  frequency  of  20  kHz.  The  beam  flux  re¬ 
flected  within  a  small  solid  angle  passed  through  a  0.633  fim  spectral  line 
filter  to  minimize  contributions  due  to  the  heating  laser  beam  reflection  and 
thermal  radiation  flux,  and  then  was  focused  onto  a  photodiode  detector. 
Simultaneous  time-resolved  surface  temperatures  on  graphite  samples  were 
measured  with  a  Si  photodetector  pyrometer  in  conjunction  with  spectral  filters 
to  minimize  reflected  1.06  fi m  laser  radiation.  Because  of  the  low  tempera¬ 
tures  at  which  Ge  and  Si  melt,  it  was  not  possible  to  measure  time-resolved 
surface  temperatures  on  these  materials. 


36Buody.  F.  P. ,  "The  P,  T  Phase  and  Reaction  Diagram  for  Elemental  Carbon, 
1979,*  J.  Geophysical  Research.  Vol.  85,  1980,  pp.  6930-6936. 

3'Ferraz,  A.,  and  N.  H.  iarch,  "Liquid  Phase  Metal-Nonmetal  Transition  in 
Carbon,"  Physical  Chemistry  of  Liquids,  Vol.  8,  1979,  pp.  289-297. 
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He/Ne  LASER 


Fig.  26.  Experimental  Apparatus 


The  surfaces  of  Ge  and  SI  samples  were  roughened  with  500  grit  paper  before 
exposure  to  insure  uniform  surface  melting  when  the  surface  was  heated  suffi¬ 
ciently.  The  pyrolytic  graphite  surfaces  roughened  to  what  appeared  to  be  a 
stable  state  as  vaporization  progressed  through  multiple  heating  cycles.  The 
Nd/glass  laser  beam  irradiated  the  material  samples  normal  to  their  surface 
and,  in  the  case  of  the  graphite,  this  heating  and  resultant  vaporization  was 
perpendicular  to  the  carbon  deposition  planes. 

Typical  oscilloscope  traces  of  the  surface  reflectance  signal  for  laser- 
heated  Ge  and  Si  are  given  in  Figs.  27  and  28.  In  these  photos,  the  upper 
trace  is  the  chopped  reflected  beam  signal,  and  the  lower  trace  shows  the  time 
history  of  the  laser  pulse  intensity.  For  Ge  and  Si  heated  to  temperatures 
below  their  melting  points,  the  reflectances  shown  in  Figs.  27a  and  27b  are 
seen  to  increase  slightly  with  increasing  temperature  and  then  to  return  to 
their  preheated  values. 

In  contrast  to  the  small  changes  in  reflectance  at  temperatures  below 
the  melting  point,  Ge  and  Si  exhibit  large  transient  changes  when  melting  does 
occur,  as  illustrated  in  Figs.  28a  and  28b.  In  these  oscilloscope  photos,  the 
time-resolved  reflectance  signal  is  again  the  upper  trace  and  the  instantane¬ 
ous  laser  heating  flux  in  arbitrary  units  is  the  lower  trace.  For  both  Ge  and 
Si,  the  reflectance  is  seen  to  increase  as  the  surface  melts  and  then  to  de¬ 
crease  to  a  value  intermediate  between  the  liquid  and  the  original  solid  sur¬ 
face  as  the  melted  phase  solidifies.  The  large  reflectance  changes  shown 
occurred  during  the  first  heating  pulse  that  visual  evidence  of  surface  melting 
was  observed  after  several  heating  cycles  with  increasing  laser  peak  heating 
flux.  The  Increase  in  reflectance  Is  due  to  the  differences  in  reflectivity 
of  the  liquid  and  solid  phases  as  predicted  by  Eqs.  (20)  and  (21)  in  Fig.  24, 
and  as  observed  during  laser  annealing  of  doped  Si  semiconductor  materials. 
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Fig.  27.  Reflectance  Below  Melting  Temperature 
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Fig.  28. 


Reflectance  at  Melting  Temperature 
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The  transient,  bidirectional  reflectance  of  roughened  pyrolytic  gra¬ 
phite  during  pulsed  laser  heating  was  over  a  wide  range  of  temperatures  by  the 
same  technique.  Figures  29a  through  29d  show  the  reflectance  response  over 
the  range  of  2850°K  to  4485*K.  These  temperatures  assume  a  surface  emissivlty 

of  0.9  at  the  pyrometer  detection  wavelength  of  0.81  /urn,  which  has  been  found 

20 

to  be  a  reasonable  value  for  graphite.  In  these  oscilloscope  photos,  the 
upper  trace  is  the  relative  reflectance  signal  and  the  lower  trace  is  the 
surface  temperature  signal.  At  the  temperatures  to  which  the  graphite  samples 
were  heated,  the  continuum  thermal  radiation  flux  from  the  heated  surface  area 
within  the  spectral  bandpass  of  the  pyrometer  becomes  greater  than  the  He/Ne 
laser  beam  flux  after  diffuse  reflection  from  the  roughened  surfaces  and  into 
the  solid  angle  viewed  by  the  pyrometer.  This  causes  an  increase  in  the  over¬ 
all  intensity  of  the  chopped  reflectance  signal  as  seen  in  Figs.  29a  and  29b. 
In  contrast,  a  change  in  surface  reflectance  is  Indicated  by  an  increased 
height  (peak  to  valley)  of  the  chopped  signal.  Thus,  it  was  necessary  to 
electronically  subtract  the  pyrometer  signal  from  the  reflected  He/Ne  beam 
detector  signal  to  obtain  adequate  sensitivity  for  the  chopped  reflectance 
signal  at  the  highest  temperatures.  This  subtraction  process  coupled  with  the 
low  level  of  He/Ne  beam  flux  from  diffuse  reflection  on  the  vaporization- 
roughened  graphite  surface  causes  the  trace  shape  shown  in  Figs.  29c  and  29d. 

As  illustrated  in  Fig.  29,  typical  pyrolytic  graphite  response  was 
consistently  observed  in  that  the  reflectance  smoothly  increased  30  to  40 
percent  during  the  time  the  surface  temperature  was  above  w3500*K,  and 
returned  to  its  initial  value  upon  cooling.  Below  maximum  surface  tempera¬ 
tures  of  «s3500°K,  no  changes  in  the  surface  reflectance  were  noted.  The 
reflectance  increase  at  the  higher  temperature  is  much  smaller  than  that  noted 
for  Si  and  Ge  when  a  melt  layer  was  created.  If  a  liquid  carbon  layer  was 
formed  on  the  roughened  surface,  it  would  be  expected  that  the  reflectance 
increase  would  be  very  large  by  virtue  of  the  reflection  surface  changing  from 
diffuse  to  specular  in  nature  and  the  inherently  greater  reflectance  of  liquid 
phases.  In  fact,  this  has  been  observed  in  reflectance  data  on  laser-heated 
roughened  iron  surfaces  that  have  not  been  included  in  this  report.  In  addi¬ 
tion  to  the  relatively  small  reflectance  change  during  heating,  post-test 
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observations  in  these  and  prevous  experiments  show  no  physical  evidence  of 
a  melt  phase  having  been  present  during  pulsed  laser  heating. 
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Fig.  29.  Reflectance  of  Pyrolytic  Graphite 
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(c)  Tma  v  =  4210°<  (  C  =  0-9) 
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(d)  Tmax  =  4485°K  {  (  =  °'9) 


Fig.  29.  Reflectance  of  Pyrolytic  Graphite  (continued) 
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From  these  considerations,  It  is  concluded  that  the  observed  small 
reflectance  increase  of  pyrolytic  graphite  probably  Is  not  due  to  melting  but, 
rather,  to  some  alternative  possibility  such  as  a  solid  phase  crystalline 
transition.  The  conclusion  that  the  carbon  triple  point  temperature  is  above 
4500 *K  is  in  contrast  to  much  of  the  earlier  experimental  results6  that 
suggested  melting  temperatures  in  the  range  of  3800  to  4200*K.  However,  as 
indicated  previously,  it  is  in  agreement  with  the  Livermore^6’^  thermo¬ 
dynamic  properties  and  thus  further  supports  our  previous  conclusions  concern¬ 
ing  carbon  vapor  pressure  and  vaporization  coefficients. 

C.  REMAINING  QUESTIONS 


The  vapor  pressure  conclusion  In  Subsection  III-B-2  is  dependent  upon 

the  assumed  emissivlty.  However,  since  the  assumed  value  of.  unity  (for  viewing 
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parallel  to  the  a-b  planes)  is  strongly  supported  by  literature  data,  ’ 
this  conclusion  should  be  considered  firm.  Likewise,  because  the  melting  con¬ 
clusion  Is  based  upon  a  direct  observational  approach,  this  conclusion  is  also 
not  likely  to  be  In  error.  The  vaporization  coefficient  conclusions  below 
4000 *K  would  also  appear  to  be  defensible  because  of  the  reliability  of  the 
emissivlty  assumption. 

For  the  pulsed  laser  data  above  4000*K,  the  pyrolytic  graphite  was 
Irradiated  and  viewed  perpendicular  to  the  a-b  planes.  The  emissivlty  In  this 
direction  for  an  ablated  surface  Is  less  well-known.  The  sensitivity  of  de¬ 
rived  vaporization  coefficients  to  the  assumed  emissivlty  value  Is  seen  in 
Fig.  17.  An  Independent  determination  of  the  pyrolytic  graphite  emissivlty 
(or  surface  temperature)  Is  needed  to  reduce  the  uncertainty  in  the  high 
temperature  vaporization  coefficient  values.  In  addition,  before  reliable 
vaporization  coefficients  for  the  polycrystalllne  graphite,  Graphnol  (see 
Appendix  B),  can  be  obtained,  It  will  be  necessary  to  determine  the  emissivlty 
of  this  carbon  as  a  function  of  temperature  for  ablated  surfaces. 
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Based  upon  the  discussion  In  Subsection  III-B-1,  relatively  strong 
dependence  of  the  pyrolytic  graphite  vaporization  coefficients  on  temperature 
may  be  the  correct  explanation  of  the  "steep"  slope  of  the  data  In  Fig.  15. 
However,  other  possible  explanations  discussed  In  Subsection  III-A-2,  and  not 
yet  eliminated,  need  to  be  given  further  consideration.  Therefore,  we  plan  to 
look  for  condensed  carbon  particles  directly  In  the  free-jet  plume  by  using  a 
He-Me  laser  and  looking  for  backscattered  radiation  from  carbon  particulates. 
Also,  we  hope  to  be  able  to  modify  an  existing  method-of-characterlstlcs  Invls 
cld  plane  computer  code  to  address  the  question  of  energy  release  effects  on 
Hach  disk  location. 
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IV.  SUMMARY  AND  CONCLUSIONS 


A  new  method,  based  upon  the  combined  use  of  carbon  laser  ablation 
data  and  analytical  modeling,  has  been  used  to  determine  carbon  vapor  pressure 
and  vaporization  coefficients  at  high  temperatures  (up  to  4600°K).  In 
addition,  a  surface  reflectivity  methou  las  been  used  to  place  a  lower  limit 
on  the  melt  temperature  of  carbon. 

The  consistently  best  agreement  of  total  vapor  pressure  data  with 
model  predictions  was  obtained  using  carbon  species  thermochemical  properties 
recommended  by  Lawrence  Livermore  Laboratory.  **  In  contrast,  comparisons  to 
the  data  with  model  predictions  using  JANAF- recommended  properties  showed  con¬ 
sistently  much  less  satisfactory  agreement.  On  the  basis  of  these  comparisons 
and  the  more  consistent  agreement  of  the  Livermore  thermodynamic  functions 
(heat  of  formation  and  entropy)  with  recent  literature  values,  It  was  con¬ 
cluded  that  the  vapor  pressure  of  carbon  Is  best  represented  by  the  Livermore 
properties.  At  4000°K,  the  JANAF- predicted  total  vapor  pressure  Is  a  factor  of 
2.5  smaller  than  the  Livermore  value.  Above  4000°K,  the  Importance  of  Includ¬ 
ing  the  contributions  of  carbon  species  Cg  and  especially  C7  to  the  total 
carbon  vapor  pressure  Is  supported  by  the  high  temperature  laser  data.  These 
species  are  Included  In  the  Livermore  publications. 

Until  now,  no  Information  has  been  available  concerning  the  tempera¬ 
ture  dependence  or  high  temperature  values  of  carbon  species  vaporization 
coefficients.  From  the  carbon  laser  ablation  data,  vaporization  coefficients 
consistent  with  the  Livermore  thermochemical  data  conclusion  above  have  been 
derived.  The  values  of  vaporization  coefficients  obtained  are  sensitive  to  the 
assumed  values  of  surface  emlsslvlty.  Thus,  particularly  for  the  data  above 
4000°K,  there  Is  a  need  to  accurately  establish  the  emlsslvlty  of  the  pyro¬ 
lytic  graphite  surface  or  to  accurately  measure  the  temperature  by  some  other 
method  (three-color  photometry).  If  it  is  assumed  that  the  high  temperature 
emlsslvlty  Is  0.97,  then  the  high  temperature  and  intermediate  temperature 
vaporization  coefficient  data  are  In  agreement  at  4000°K. 
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An  auxiliary  He-Ne  laser  was  used  to  monitor  the  carbon  surface  re¬ 
flectivity  as  a  function  of  time  during  the  laser  pulse.  NASA  data  at  measured 
temperatures  up  to  4500*K  showed  no  evidence  of  the  large  Increase  in  reflec¬ 
tivity  expected  when  the  surface  melts.  Such  changes  (Increases)  were  success¬ 
fully  observed  for  other  group  IV-A  elements.  Thus,  it  was  concluded  that  the 
melt  temperature  of  carbon  Is  greater  than  4500*K. 

The  conclusions  regarding  vapor  pressure  and  vaporization  coefficients 
are  Important  because  the  mass  loss  rate  Is  proportional  to  the  Langmuir  pres¬ 
sure,  l.e.,  the  product  of  the  vaporization  coefficient  and  the  vapor  pressure. 
In  addition,  the  melt  temperature  conclusion  Is  technologically  important, 
since  melting  could  reduce  the  energy-dissipating  efficiency  by  as  much  as  a 
factor  of  three  If  the  liquid  carbon  Is  stripped  away  without  vaporizing. 
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APPENDIX  A 


DATA  SUMMARY 


The  experimental  data  analyzed  in  this  report  are  stmmarized  in 
Tables  A-l  and  A-2  ( CVI  laser  data)  and  in  Figs.  A-l  through  A-4  (pulsed  laser 
data).  A  more  complete  description  is  given  in  Refs.  10  and  38.  The  data 
lines  in  Fig.  A-4  were  obtained  by  combining  the  information  in  Figs.  A-l 
through  A-3. 


3®Cov1ngton,  A.,  "Ames  Free-jet  Laser  Vaporization  Experiments,' 
NASA  Ames  Research  Center,  Moffett  Field,  CA,  September  1980. 
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Table  A-l . 


Pyrolytic  Graphite  Ablation  Data* 


A  H 

8  I  £  H  T 

P  *  E  S  S  U 

R  £ 

Laser  Power 
<k«) 

Torr 

Atm 

0.5 

1 

0.005 

0.02 

0.05 

0.1 

0.2 

0.25 

0.5 

1.0 

1.59 

3384V 

1.  £61119/ sec 

1.89 

34  5  5V 

2. 5 lug /sec 

2.18 

3687  V 
17.2«g/sec 

2.32 

2.47 

3560V 

7.80»g/sec 

2.83 

3592V 

9.41ng/sec 

3765V 

38.6mg/sec 

2.90 

3810  V 

21 .7ng/sec 

3.05 

3587V 

8.67nj/$ec 

3.20 

3642  V 
13.3«g/sec 

3.34 

3631 V 
13.9mg/sec 

3687  V 
20.2i»g/sec 

3698  V 
18.8wg/sec 

3740V 

25.0i*g/sec 

3790V 

29.4«g/sec 

3830V 

31. 2»g/ sec 

3880V 

34. 1  mg/ sec 

3.34 

3731V 

18.7«g/sec 

3.34 

3720V 

17. 2*9/ sec 

3.79 

3760V 

37.fc*g/sec 

3.92 

3912V 

47.9ng/sec 

4.10 

3810V 

47. lug/ sec 

3822  V 

68. 3ag/ sec 

3878V 

59.4«g/$ec 

4.58 

3884  V 
72.9mg/sec 

_ 

♦The  nunbers  given  In  each  column  are  the  measured  surface  temperature  and  mass 
loss  rate. 


Table  A-2.  Graphnol  Ablation  Data* 


User  Power 

U») 

A  M  8  I  E 

NT  P  R  E 

S  S  U  R  E  ) 

Torr 

_ _ 1 

1 

0.05 

0.1 

0.2 

0.5 

1.0 

1.05 

1.45 

2870V 
1.11*9/ see 

1.96 

3335V 

1 .89»g/sec 

2.04 

3217  V 

5. 11*g/sec 

2.10 

3400V 
18.0*g/ sec 

2.33 

3466  V 

4. 90»g/sec 

2.41 

3521V 

18. 2*g/ sec 

2.76 

3532V 

18.0ag/sec 

2.83 

3535V 

12.5ag/sec 

2.83 

3530V 

9.33ag/sec 

2.90 

3565V 

18. lag/ sec 

3715  V 

30. Bag/ sec 

2.91 

3S32*K 

8.52ng/sec 

2.98 

3565V 

11 .7ng/sec 

2.98 

3526V 

10. lag/sec 

2.98 

3687V 

30.6ag/sec 

3.05 

3565V 

12.6ag/sec 

3.05 

3520V 

11.3ag/sec 

3.05 

3570V 

21.4ag/jec 

3.05 

3648V 

27.3ag/sec 

♦The  numbers  given  In  each  column  are  the  measured  surface  temperature  and  mass 
loss  rate. 
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Table  A-2.  Graphnol  Ablation  Data  (Continued) 


unr  Power 
<»») 

A  M  8  I  E  N 

T  PRESSURE 

Torr 

At* 

1 

0.05 

0.1 

0.2 

0.5 

1.0 

1.05 

3.19 

3681*8 

29.4ag/sec 

3.34 

3609*8 

27.0ag/sec 

3.41 

3631*8 

29. 5»g/ sec 

3.63 

3754*K 

37.2*g/tec 

3.99 

3777*8 

42.9ag/sec 

4.07 

3600*8 

31 .6*9/ sec 

4.21 

3794*8 

49. 4*0/ sec 

4.29 

3587*8 

31.2ag/sec 

4.36 

3687*8 

25.2ag/sec 

4.36 

3665*8 

46.8ag/sec 

4.36 

3642*8 

39.3ag/sec 

4.36 

3687*8 

42,4ag/sec 

4.36 

3833*8 

56.5ag/sec 

4.50 

3743*8 

68. Bag/ sec 

4.58 

3830*8 

52.6ag/sec 

4.58 

3833*8 

47 .  Sag/ sec 

5.30 

_ 

3884*8 

1 1 7ag/sec 

P  (atm) 


(W/cm  ) 


ASSUMED  €  ■  0.9 


EXP 

GRAPHITE  DATE 


•  6  PYROLYTIC  (super  temp)  1980 

▲  ATJ  1980 

+  GRAPHNOL  1980 

■  ATJ  1976 


Fig.  A-3.  Peak  Surface  Temperature  versus  Peak  Laser 
Flux  for  Vaporization  In  Argon 


APPENDIX  B 


GRAPHNOL  DATA  INTERPRETATION 


As  seen  In  the  data  summary  given  in  Appendix  A,  a  large  quantity  of 
experimental  data  were  also  obtained  for  the  fine-grain  polycrystalline 
graphite,  Graphnol.  Because  the  Graphnol  thermal  conductivity  is  nearly  iso¬ 
tropic,  the  Graphnol  cylindrical  samples  were  heated  to  a  bright  glowing  orange 
by  axial  conduction  away  from  the  laser  Irradiated  area.  This  caused  mass  to 
be  lost  from  the  cylindrical  surface  due  to  heterogeneous  reaction  with  oxygen 
In  the  ambient  chamber  atmosphere.  A  considerable  effort  was  devoted  to  elimi¬ 
nating  the  oxidation  mass  loss  contribution  to  the  total  mass  loss  so  that  the 
sublimation  mass  loss  rates  could  be  used  to  obtain  vapor  pressure  and  vapori¬ 
zation  coefficient  Information.  Unfortunately,  no  satisfactory  unambiguous 
method  of  doing  this  was  found,  even  for  the  very  low  pressure  (1  Torr)  data. 

The  1  Torr  data  points  from  Appendix  A  are  plotted  In  Fig.  B-l.  Also 
shown  is  the  pyrolytic  data  line  from  Fig.  10.  Note  that  both  lines  have  about 
the  same  slope.  Interpreting  these  data  directly,  as  discussed  In  connection 
with  Figs.  10  through  12  in  the  report,  we  conclude  that  the  vaporization  coef¬ 
ficient  for  Graphnol  is  greater  at  a  given  temperature  than  that  for  pyrolytic 
graphite  by  a  factor  of  about  1.7.  The  amount  of  mass  lost  by  oxidation  Is  not 
known,  and  the  emlsslvlty  of  the  Graphnol  is  likely  to  be  less  than  unity  (as 
assumed  in  the  data  reduction).  Since  correction  for  both  these  effects  would 
reduce  the  vaporization  coefficients  derived  from  the  data.  It  is  probable  that 
the  vaporization  coefficient  for  Graphnol  is  actually  less  than  that  for  pyro¬ 
lytic  graphite  at  a  given  temperature. 
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The  high  temperature  Graphnol  data  from  Appendix  A,  when  Interpreted 
as  discussed  in  the  report  (see  discussion  of  Figs.  13  through  16),  give 
vaporization  coefficients  as  seen  In  Fig.  B-2.  Again,  shown  for  comparison  Is 
the  pyrolytic  graphite  data  line  from  Fig.  16.  If  the  emlsslvity  of  Graphnol 
and  pyrolytic  graphite  are  about  the  same,  then  these  results  Indicate  that 
carbon  vaporization  coefficients  for  Graphnol  are  approximately  a  factor  of  5 
smaller  than  those  for  pyrolytic  graphite  at  a  given  temperature. 


Fig.  B-l .  Comparison  of  Low  Pressure  Pyrolytic  Graphite 
and  Graphnol  Hass  Loss  Rate  Data 


APPENDIX  C 


INFORMATION  AVAILABLE  FROM  CW  VERSUS  PULSED  LASER  DATA 


As  Indicated  schematically  in  Fig.  8,  carbon  vapor  pressure  informa¬ 
tion  is  obtained  only  from  the  CW  laser  data.  The  purpose  of  this  Appendix  is 
to  discuss  why  vapor  pressure  information  cannot  be  derived  from  the  pulsed 
laser  data. 

The  theoretical  model  predictions  (Subsection  II-A)  can  best  be  sum¬ 
marized  by  referring  to  Fig.  20.  When  the  ambient  pressure  pa  is  much  less 
than  the  carbon-saturated  vapor  pressure  p$  at  temperature  T  ,  then  the 
predicted  mass  loss  rate-surface  temperature  behavior  is  described  by  the 
straight  line  in  Fig.  20a.  From  the  CW  laser  data  shown  in  Fig.  20a  or  the 
pulsed  laser  data  shown  in  Fig.  15,  the  Langmuir  vapor  pressure  can  be  obtained 
via  Eq.  (11).  Then,  once  p$  is  obtained  by  independent  means,  Eq.  (11)  can 
be  solved  for  the  vaporization  coefficient  a. 

When  the  above  condition,  l.e.,  pa«pe,  is  not  satisfied,  the  pre- 
dieted  mass  loss  rate-surface  temperature  behavior  is  described  by  the  various 
curved  lines  (each  for  a  different  pm)  shown  in  Figs.  20b  and  20c.  In  Figs.  21 
through  23,  CW  laser  data  are  compared  with  model  predictions  using  different 
thermochemlcal  data  sets  to  determine  the  "more  nearly  correct"  vapor  pressure 
relationship.  In  order  to  do  this,  it  is  necessary  that  the  data  be  In  the 
curved  lines  region  (see  Fig.  20b)  of  the  m  -  1/T  plane.  Since  the  pulsed 
laser  data  (Fig.  15)  all  lie  In  the  straight  line  region  of  model  predictions, 
the  method  applied  to  the  CW  data  In  Figs.  21  through  23  to  obtain  vapor  pres¬ 
sure  information  cannot  be  used. 


The  reason  that  the  pulsed  laser  data  all  lie  along  a  straight  line  is 
fundamentally  inherent  to  the  experimental  method,  i.e.,  stagnation  pressure  pQ 
is  derived  from  measured  Mach  disk  geometry.  To  have  a  Mach  disk,  the  flow 
must  be  sonic  (choked)  at  the  surface  and  expand  in  a  supersonic  free- jet 
expansion.  The  condition  required  for  sonic  flow  at  the  surface  is  pa«  p$. 
From  above,  recall  that  this  is  the  same  condition  as  that  required  for  the 
straight  line  m  versus  1/T$  behavior.  Thus,  all  the  pulsed  laser  data  are  in 
the  straight-line  region. 
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